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Renormalization of nonequilibrium dynamics in FRW cosmology 


Jurgen Baacke *, Katrin Heitmann 1, and Carsten Patzold 1 
Institut filr Physik, Universitat Dortmund 
D - ^^221 Dortmund, Germany 
(June 10, 1997) 

We derive the renormalized nonequilibrium equations of motion for a scalar field and its quantum 
back reaction in a conformally flat Friedmann-Robertson-Walker universe. We use a fully covariant 
formalism proposed by us recently for handling numerically and analytically nonequilibrium dynam¬ 
ics in one-loop approximation. The system is assumed to be in a conformal vacuum state initially. 
We use dimensional regularization; we find that the counter terms can be chosen independent of the 
initial conditions though the divergent leading order graphs do depend on them. 
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I. INTRODUCTION 

Nonequilibrium processes in cosmology have been considered recently by various authors. The main interest has 
been centered around the possible inflationary period of the universe (see e.g. jl| ||]) and the subsequent reheating 
H,|| • It has already been found by considering the parametric resonance @ associated with oscillations of the inflaton 
field and by exact computations including back reaction both in Minkowski space [Q and in an expanding universe 
||[j) that the time-dependent inflaton field produces particles or classical fluctuations preferentially at low momenta 
and not in a distribution corresponding to thermal equilibrium. This process of particle production has therefore been 
termed preheating fUTof. 

The equations of motion for nonequilibrium systems have been presented by various authors |TT] |l3j using the CTP 
formalism introduced by Schwinger ||l4|| and Keldysh 0- Their application to inflation within a conformally flat 
FRW universe has been initiated by Ringwald [ 16] [jl|; they have been recently implemented numerically in jl9|] and 
p(i| . Similar computations have been performed recently in configuration space for the case that the fluctuations are 
treated as classical ones pllp2]]. Again fluctuations with rather low momenta are strongly excited, thus justifying the 
classical approximation. Wpart form such exact numerical computations there exist also various analyses based on 
analytical approximations to the solution of the mode equations |)|,^3] ^7| . 

We have proposed recently a computation scheme for nonequilibrium dynamics which has two aspects: on the 
one hand it separates cleanly the divergent and the finite contributions of the quantum fluctuations in one-loop 
approximation, and so is attractive for numerical computations. On the other hand it leads, by the analysis of the 
divergent leading order contributions, to a simple formulation of the renormalized equations of motion. The fact that 
the divergent contributions are removed from the numerical computation allows a free choice of regularization. We 
have chosen dimensional regularization since it is easier to handle in the presence of quartic and quadratic divergencies 
than e.g. Pauli-Villars regularization.This scheme has been applied and implemented numerically for scalar fields | p8[ , 
fermion fields and for the SU(2) Higgs model |2l|. We will present here the renormalized equations of motion for a 
scalar field in flat Friedmann-Robertson-Walker (FRW) cosmology. 

We will not attempt here a numerical implementation. There exist already several numerical analyses based on 
various assumptions on the initial conditions and on the parameters of the theory. We do not expect that our 
formulation will lead to major differences in the qualitative features of the results. Nevertheless, we think that in the 
presence of quartic and quadratic divergencies it is important to formulate the renormalization scheme in a way that 
is fully covariant and independent of the initial conditions. It turns out that the resulting formalism represents at the 
same time a rather attractive numerical computation scheme. 
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II. FRW COSMOLOGY 


We consider the Friedmann-Robertson-Walker metric with curvature parameter k = 0, i. e. a spatially isotropic 
and flat space-time. The line-element is given in this case by 

ds 2 = dt 2 — a 2 (t)dx 2 . (1) 

The time evolution of the a(t) is governed by Einstein’s field equation 

+ aHj^J + + Ag^ u = —k(T m1/ ) , 

with At = 8irG. The Einstein curvature tensor G^ is given by 

Gin/ — R^iu 2 9^R • 

The Ricci tensor and the Ricci scalar are defined as 

R[lv = R/j. v \ i (4) 

R = tf v R» v , (5) 

where 


( 2 ) 

( 3 ) 


R 


\ 


= - d a T* a T a a0 


pA p(7 
-*■ acr 1 - 7^0 • 


( 6 ) 


The terms H$ and H^J arise if terms proportional to R 2 and R^R^ are included into the Hilbert-Einstein action. 
If space-time is conformally flat, these terms are related by 


if( 2 ) = i H (1) 
±J - fiu 3 fiu 


( 7 ) 


so that we can set f3 = 0 in (|j) without loss of generality f]f| . We also replace H$ by in the following. 

These terms are usually not considered in standard cosmology. They are included here, as well as the cosmological 
constant, only for the purpose of renormalization; they will absorb divergencies of the energy momentum tensor. So 
in principle they should appear on the right hand side as counter terms for the energy momentum tensor; they are 
related to the coefficients of these counter terms [] by A = k5A and a = nSa. We will need a further counter term 
proportional to G which can be considered as a wave function renormalization for the gravitational field; either we 
have to replace G by ZG on the left hand side of the Einstein equations or to introduce a counter term SZG„„ 
for the energy momentum tensor. Again both alternatives are related by Z = 1 + n5Z. Writing the renormalized 
equation (§ in the form 

(jl+sz) G a- + (“ + te ) H a- + (a + SA) g^ = T (8) 

we can identify SZ as a correction to the gravitational coupling via At -1 — > At -1 + SZ. As usual, we can reduce the 
Einstein field equations to an equation for the time-time component and one for the trace of G^j,, the Friedmann 
equations 


G t t + oiHtt + A — —At Ttt , 
G£+aff£+4A=-/tT£. 


For the line element (|1,) the various terms take the form [|30| 


Gtt(t) 

G*{t) 

H u (t) 

H£(t) 


-3 H 2 (t) 
-R(t) , 


-6 

-6 


H(t)R(t) + H 2 {t)R(t ) - -R 2 (t) 


R(t) + 3 H(t)R(t) 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


'See 


for their precise definitions. 
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with the curvature scalar 


and the Hubble expansion rate 


Rpt) = 6 H{t) +2H 2 (t) 


H(t) = 


«(*) 

a{t) 


(15) 


(16) 


III. NONEQUILIBRIUM EQUATIONS FOR THE SCALAR FIELD 

The Lagrangian density of a c(> 4 -theory in curved space-time is given by 

C = - im 2 $ 2 - §R,$ 2 - fp$ 4 } , (17) 

where R{x ) is the curvature scalar and f the bare dimensionless parameter describing the coupling of the bare scalar 
field to the gravitational background. We split the field $ into its expectation value </> and the quantum fluctuations 
ip: 


with 


$(f,t) = (p{t) +ip{x,t) , 


(pit) = (®{x,t)) = 


Tr<hp(t) 
Tr pit) 


(18) 


(19) 


where pit) is the density matrix of the system which satisfies the Liouville equation 

i^- = m),pit)}. ( 20 ) 

The one-loop equation of motion of a scalar field with A <p 4 interaction has been obtained in the FRW universe by 
Ringwald [Q; we follow closely his formulation. The equation of motion for the classical field is 

<p + 3 Hj> + (m 2 + £R)<P + Kb 3 + ^-(ip 2 )& = 0 . (21) 

6 2 

The expectation value of the quantum fluctuations pip 2 ) can be expressed as 

pip 2 ) = — iGpt, x; t, x) (22) 


in terms of the non-equilibrium Green function G(t, x\ t', x r ) which satisfies 


-g^2 + + a + m2 + £-^(i) + 


Git,x;t',x') = - 37 -t- 5it,x;t',x') . 
a 6 it) 


(23) 


The boundary conditions for this Green functions will be given below. Due to the presence of the term Hit)d/dt 
the differential operator on the left hand side of this equation is non-hermitian. It is made hermitian by introducing 
conformal time and appropriate scale factors. Conformal time is defined as 


t = ■ <24) 
o 

In conformal time the line-element (jl|) reads 

ds 2 = G(r)(dr 2 — dx 2 ) , (25) 

where the conformal factor Cpr) is given by 
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C{t) = a 2 (r) . (26) 

We further rescale the scalar field and its quantum fluctuations by introducing the dimensionless ‘conformal’ fields 


(p(r) = a(t)<j)(t ) , 

= a(t)ip(t,x) . 

The Green function is rescaled accordingly via 

G(x,t; x 1 ,t') = a(t)a(t')G(x, t; x, t') . 

The equation of motion of the classical field < p(r) now becomes 

ip" + a 2 [m 2 + (£ - l) -R] ip + + ^<pG(r, x; r, x) = 0 , 

where the primes denote derivatives with respect to conformal time. The two-point-function G now satisfies 

d 2 


dr 2 


- V 2 + M 2 (r) 


G{x,t\ x,t') = —5(x,t',x',t') . 


Here M 2 (t) denotes the square of the effective mass term the fluctuation field %fi(r,x) 


M 2 (t) = C{t) 


m 


+ (£ ~ g) R ( t ) + o 


\<P 2 (t) 


2a 2 (r)J ' 


(27) 

(28) 


(29) 


(30) 


(31) 


(32) 


The problem of determining the Green function is now essentially reduced to the equivalent problem in Minkowski 
space. We expand the fluctuation field in terms of the mode functions U k {r) exp (ikx) via 


4>(t,x) = 


' d 3 k 


c(k)Uk(r) exp (ikx) + c 1 (k)U k (r) exp {—ikx) 


The functions t4(r) satisfy the mode equation 

C/0(r)+H2( T )C/ fe (r)=O, 

with 

^l(r) =k 2 + a 2 (r) [m 2 + (£-£) R(t)] + 2 (r) . 


We further impose the initial conditions 


with 


C4(0) = 1 ■; £4(0) = -ififc(O) , 


fifc(O) = \A 2 + M 2 (0) . 


(33) 

(34) 

(35) 

(36) 

(37) 


In the following we will use the short notation = 12^(0). The nonequilibrium Green function Gi c (t,x;t',x') can 
by expressed in terms of the mode functions via 


G k (r,x; t', x') = 


dk 1 - T')U k {T)U k {T')exp{ik{x - x')) 

+9{t’ — r)U k (r f )U k (r) exp(— ik(x — f'))} • 


(2tt) 3 2f2fc 0 


The expectation value of the fluctuation fields is given, therefore, by the fluctuation integral 

' d 3 k |t/fc(r) | 2 


jC(r) = (ip (r)) = —iG(r, x\ r, x) = 


(2t r) 3 212 


kO 


(38) 


(39) 
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The unrenormalized equation of motion of the inflaton field reads 

<P" + a 2 0) [m 2 + (£ - i) R(r)] ip + + ^pR(r) = 0 . (40) 

The regularization of the fluctuation integral and the renormalized form of this equation will be discussed below. 

The boundary conditions for the mode functions and the corresponding definition of the fluctuation integral are 
related to the definition of the initial quantum state and/or density matrix of our nonequilibrium expansion. It would, 
e.g., look different if we had used the time variable t instead of conformal time r. These initial conditions, a conformal 
vacuum state |j(J , have already been chosen in jlTj , they are discussed to some extent in ||. The initial conditions for 
the cosmological expansion are even more subtle since even in principle we cannot start from an equilibrium situation, 
except if the classical field is at a minimum of the effective potential and the energy momentum tensor vanishes there 
- obviously not an interesting situation. When considering nonequilibrium processes in Minkowski space-time one can 
imagine that an initial classical field is maintained away from the minimum of the effective potential by a source (like 
capacitor plates for an electric field) which is switched off at time zero. In such a situation it makes sense to assume 
that the initial state is the vacuum state corresponding to that initial classical field. In the cosmological context the 
presence of a nonvanishing classical energy leads already to a cosmological expansion which in turn influences the 
quantum fluctuations. So starting with the conformal vacuum state or a thermal state constructed on it seems a priori 
artificial. Since there is no solution of principle to this problem one has to hope (or try out by numerical experiments) 
that the precise initial conditions do not influence too strongly that period of the cosmological expansion one wants 
to study. 


IV. THE ENERGY MOMENTUM TENSOR 


In order to formulate Einstein’s field equation we have to discuss the energy momentum tensor of the scalar field 
in curved space time. For a classical field it reads Q 


T„u = (1 - 2 + ( 2 £ - \) g^g pa (!).p(j). a - 2 

+2 £g,jv<f>n<f> - £,Gp U (j) 2 + |m 2 s F i) 2 + . 


(41) 


In the conformally flat FRW metric the energy momentum tensor is diagonal. One obtains for its time-time component 
and its trace 


T t c i = - £G tt (j) s + 6 iHU 


it = X 2 + + lx 4 

Tp cl = -X + 2m 2 X + + 6£(# + 0 2 + 3 H<M>) . 


(42) 


We introduce again conformal time and the conformal rescaling of the fields. Furthermore, we include the quantum 
fluctuations of the field ip. We obtain 


1 


1 


A 


<P + + (! - 6 0 1 7T3 <P - —W' 


2 a 4 

1 


2 a 2 
d 3 k 1 


H 2 


H 


2a 2 ' 


i 2 J (27r) 3 2fl 2 0 | 2a 2 


\u' k \ 2 - 


2 a 2 


n(T)\u k \ 2 -±(z-l)R\u k \ 2 


H d 


-hm - m 2 \u k \ 2 + m - i)--h^ i 


a dr 


( 43 ) 


and 


2 We continue to consider Ttt instead of T tt = a 2 Ttt for convenience. 
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, J 2 tt2 

ip H 2 


77 


T n = (! - 60-r- T p + 2—<p<// + 6-7^" +2m —+ 


£ 


1 /• d 3 fc 1 


a 2 J (27r) 3 2S2 


fcO 


(1 - 60 




2 <A 2 , A 4 

-m 


m 2 


2 + ^""( r ) 2 
ar a z 


, 2 ra 2 


— + imi 2 --/W 

ar a ar 


A <y 9 2 




i%i ! 


Energy density and pressure are related to the energy momentum tensor via 


( 44 ) 


T tt = £, 

Tjl = £ - 3p . (45) 

It is straightforward to show, using the equations of motion for the classical field and for the mode functions (|2ij|), 
that the energy is covariantly conserved: 

£'{r)/a{r) + 377(r)(p(r) + £ (r)) = 0 . (46) 


V. PERTURBATIVE EXPANSION 


In order to prepare the renormalized version of the equations given in the previous section we introduce a suitable 
expansion of the mode functions, which was successfully used in p^,^9[ for the inflaton field coupled to itself, to 
fermions, and gauge bosons in Minkowski-space. Adding the term M^JO) on both sides of the mode function equation 
it takes the form 


dr 2 


^fcO 


Uk{r) = -U(r)C/ fc (r) , 


(47) 


with 


V(r) 


M 2 {t) — M 2 (0) , 

r _ -i 1/2 

k 2 + M 2 (0) 


(48) 


(for the definition of M 2 (r) see eq.(|32|)). The first and second derivatives of this potential will be needed in the 
perturbative expansion and are given by 


V'{t) = 2clHM 2 (t) + a 2 (£ — A) R' + \<p<p' — XaHtp 2 , 
V" = \a 2 RM 2 - 2a 2 H 2 M 2 + 2aHV'( t) 

+2a 3 H (£ - A) R' - 2XaHpxp' + a 2 (£ - A) R" 

+X(ip<p" + ip' 2 ) — X—Rp 2 + Xa 2 H 2 ip 2 . 
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Including the initial conditions (36) the mode functions satisfy the equivalent integral equation 


(49) 


(50) 


t4(r) = e~ in *° T + j dr'A fc , ret (r - r')V(r')U k (T') , (51) 

0 


with 

Afc !re t(T - r') = -2— 0(r - T / )sin(f7 fc o(r - r')) . (52) 

We separate Uk{r) into the trivial part corresponding to the case V(t) =0 and a function hk(r) which represents the 
reaction to the potential by making the ansatz 


6 















E4(r) = e~ in ^ T (l + h k (r)) . 


(53) 


hk(r) satisfies then the integral equation 

r 

/ife(r) = Jdr'Ak,ret (t - + /i fe (r , ))e in,cCl(T_r } , 


(54) 


and an equivalent differential equation 

K( t ) - 2 iflkoh' k {T) = -V(r)(l + /i*(r)) , (55) 

with the initial conditions hk{ 0) = hk( 0) = 0. We expand now /ifc(r) with respect to orders in V(t) by writing 

hk (t) = /if } (r) + /if } (r) + /if 3 (r) + • ■ ■ (56) 

= h k\ T ) + h k\ T ) ^ (57) 

where /if E) is of n’th order in V{r) and /if 3 ("r) is the sum over all orders beginning with the n’th one: 


cv)=emtv). 


(58) 


The /if 3 are obtained by iterating the integral equation (|54|) or the differential equation (|55|). The function /if 3 (r) 

is identical to the function hk (r) itself which is obtained by solving ( p5| ). The function /if 3 (r) can again be obtained 
by iteration via 


/if (r) = J dr' A fc , ret (r - T0H(r')/if V)^ 0 ^' 3 , (59) 

0 

or, using the differential equation, via 

h "k ] ( T ) - 2if} fc0 /i'f 3 (r) = ~V (r)/if 3 (r) . (60) 

This iteration has the numerical aspect that it avoids computing /if 3 via the small difference /if" 1 — /if 3 . However, the 
integral equations are used as well in order to derive the asymptotic behaviour as flko —> oo and to separate divergent 
and finite contributions. We will give here the relevant leading terms for /if 3 (r) and We have 

r 

h { k\r) = J dr , (exp(2/H fe0 (r - r')) - l)H(r') . (61) 

0 


Integrating by parts we obtain 


h£\r) = - 1 f - 1 1 




c/r H(r') - F(r) + —3- / dr' exp(2/H fc0 (r - t'))F'(t') , 


fcO 


‘'fcO 


or, by another integration by parts, 


h { k\r) = 


2n 


kO 


VF(r,) - A v(t)+ w/' {T) 


lo 


dr' exp(2iflfc 0 ( T — T ')W"( T ') ■ 


(62) 


(63) 

(64) 
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We will need often the real part of li k 1 for which we find 

r 

R e /i£ 1) (r) = --^r V( y T ) + 4^2- Jdr' cos{ 2 VL k0 {T - t'))V\t') 


4 ft 


^ y(r) + 4 ^ C(y, ’ T) ' 


fco 


( 65 ) 


feO 


Since various Fourier integrals will appear in the finite parts of the fluctuation integrals we introduce for later conve¬ 
nience the notations 


T 

«5(/,r) = J dr / sin(2fl i: o(T-r / ))/(r') , 
0 

T 

C(/,t) = J dr'cos(2f2fc 0 (r — t , ))/(t / ) . 


( 66 ) 

(67) 


For the leading behaviour of h k (r) we find 


r c 

4 2) M = -774r /W dr"V(r')V(r") + O(il k0 3 ) . 


( 68 ) 


o o 


In terms of this perturbative expansion we can write the mode functions appearing in the fluctuation integrals in the 
equation of motion and in the energy momentum tensor as 


|C4| 2 = l + 2Re4 1) + l4 1) | 2 - 


(69) 


and 


\U k \ 2 — (l + 2Re/4 ^ + | h k + | h k ^| 2 

—iClko ^—2 ilmh'J^ — 2ilmh^*h'^'j . 

As the potential is real, the leading behaviour of the sums is 

1 + 2Re/4 T) + \h®\ 2 = 1 - ^tV(t) + ^ sin(2fI feo r)F'(0) + ^R"(r) 

/cO kO /cO 

--^(20 fco r)I/"(0) + -^(r) + O(Q k0 5 ) , 


(70) 


(71) 


L k0 


and 


- 2iImh'P - 2ilmh ( i ) *h'i 1) = -- —L- sin(2f2 fco r)y'(0) - -i- V"(r) 




fco 


2Q 2 




3i 


From the Wronskian relation 


we obtain the relation 


2if2fco ^ 


+ —^ cos(2O fc0 r)y"(0) - T 7 ^y 2 (r) + O(Q k0 4 ) 
^ L k0 ^ L k0 

u k u* k ' - u' k u* k = 2in k0 

2Re4 1} + l4 1} | 2 ) ^ 2ilm/if } - 2iIm/4 1) X (1) = 0 , 


(72) 

(73) 

(74) 


which proves to be useful in simplifying the mode integrals occuring in the energy momentum tensor. 



VI. RENORMALIZATION OF THE EQUATION OF MOTION 


The expansion of the mode functions allows a free choice of the regularization scheme. We will use dimensional 
regularization here since in the presence of quadratic and quartic divergencies it is much easier to handle than in 
Pauli-Villars regularization. We will show that in this scheme the divergence structure of the equation of motion and 
of the energy momentum tensor have the correct form for a consistent renormalization. Furthermore, we will present 
the explicit form of the finite parts of the fluctuation integrals occuring in the equation of motion and the energy 
momentum tensor. 

The fluctuation integral of the equation of motion can be split into a divergent and a convergent part. Using (|65| ) 
we obtain 


Hr) - / 


d 3 k 1 + 2Reh^(r) + |h^(r)| 2 

WF 20 ^ 

’ d 3 k 1 (. 1 

(2t r) 3 2kl k0 


1 - 


V(r) 


r C(V\t) 


2f2feo 2f2fco 

+2Re4 2) (r) + |4 1} (r)| 2 ) . 


(75) 


The first two terms in the integrand have to be regularized. We first rewrite the basic equation of motion, including 
appropriate counter terms, as 


<p" + a 2 [m 2 + 5m + (£ - | + <5£) R]<p+ —V^ 3 + ~ 0 ■ 

Next we separate from the term \(pT / 2 the dimensionally regularized divergent parts 


<P(r) 


d 3 k 1 


(2ir) 3 2Q fc0 


= P e ^P( T ) J 


d d k 

(2'k)' 1 


k 2 + M 2 { 0) 


1/2 


AM 2 (0)ip(r) / 2 47r/i 2 

e Af 2 ( 0) 7 


32t r 2 


(76) 


(77) 


and 


- jw ) 5 W7 },.« = j 


d d k 


( 27r ) d \p + M 2 (0)j 3/2 


A <p(t)V(t) f 2 47r/i 2 

7 +ln m 0)- 7 


32tt 2 


(78) 


Recalling that V(t) = M 2 (r) — M 2 ( 0) and M 2 (r) = a 2 (r) [m 2 + (£ — l/6)ii] + \<p 2 /2 one sees that M 2 (0) cancels 
for the divergent terms and that therefore the counter terms can be chosen independent of the initial conditions as 


Sm 2 = 

A m 2 J 

f 2 

+ In 

47 Tfl 2 
- 0 -7 


32tt 2 1 

u 

m z 

S\ = 

3A 2 j 

f 2 

+ In 

47r/i 2 
- 0 -7 

32t r 2 | 

u 

ra z 


1 

II 

1) 

/ - + In 4?r/i2 

32tt 2 

l 6 

m z 


(79) 

(80) 
(81) 


The renormalized equation of motion now reads 


<p" + a 2 


to 2 + A to + (£ + A£ - -)R 


A T AA o A 

V? H-g —P + 2 <PSfin = 0 , 


(82) 


with 
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(83) 


fin — 


M 2 ( 0) 

I 67 r 2 


d 3 fc 1 


and with the finite corrections 


/ (27t) 3 2f2fc 0 V2f2 fe0 

a 2 ^ Tr)2 

A 777 2 
A A 

a? = 


I C(V',T) + 2Reh® + \h£ ) \ 2 


32tt 2 

n M 2 (0) ’ 

3A 2 

m 2 

32tt 2 

n M 2 (0) ’ 

A(€- 

■*) h 


32t r 2 


M 2 (0) 


(84) 

(85) 

( 86 ) 


VII. RENORMALIZATION OF THE ENERGY MOMENTUM TENSOR 


In order to derive a renormalized form of the Friedmann equations we have to renormalize the energy momentum 
tensor as well. In principle this has been discussed long ago and it will not be a surprise that the divergent parts are 
in one-to-one correspondence to those given e.g. in |30j. However, we have to discuss this subject in the framework of 
nonequilibrium quantum field theory, and we are interested in particular in the precise form of the finite parts which 
will be the subject of a numerical computation. In order to renormalize the energy we introduce the available counter 
terms into the unrenormalized expression so that it reads now 


£ = —7 ip r 

2o 4 Y 


1 

2a 2 


(m 2 


Sm 2 )ip 2 

T2 


A + S A . 


-6 (£ - i + O + <5A + SaH tt + 5ZG 

/A >—... . n 1 1 . ( 1)|2 


■ d 3 k 

(^j 3 


J2( 1 +2R eA ® 


I'*?! 2 ) 


+ 


l^fl 2 


2a 2 V k 1 k 1 ) 4fi fc0 a 2 

-4^2 (~2*Im h’^ - 2Hm/4 1) X (1) ) 

+ 3o^ y<T) ( 1 + 2Re '‘‘ Cr, + A m l 2 ) 
1 


4H fc0 (6 ^ 1} V 6 


1 


4H 


fcO 


i 2 ^j (l + 2Reh® + \h®\ 2 ) 

(6^ - !) —— (i + 2Re/i^, 1) + J/i^l 2 ) } . 


(87) 


The divergent parts of the fluctuation integral are 

1 f d 3 k 


£ 


div, flue — 9 

a 2 


(27r) 3 


— +- : - V(t) - 

2a 2 T An ~2 y T ) 


1 


4f2fcoa 


16H 3 


;V 2 (t) 


kO u 


-3(6f-l)( e + ff2 


R 


' 2Hfco 




V(r)) 


-i(6£-i) 


Dimensional regularisation of the divergent integrals yields 


H 1 

a 4H^ 0 


V\t) 


( 88 ) 


d 3 fc 1 


n 


— + -A-U(r) - 

o in ^ ' 


1 


(27r) 3 a 2 2 
M 4 (t) f 2 
647r 2 a 2 


4Hfco 

- + i n iE^_ 

e M 2 (0) 


16H^ 0 


U 2 (r) 


M 4 (0) M 2 (0)M 2 (t) 


[m 2 + (ij~\)R+^ 2 Y 

647T 2 


1287r 2 a 2 

2 47 t ft 2 

e n M2(0) 


327r 2 a 2 


( 89 ) 


- 77 + 


M 4 (0) M 2 (0)M 2 (r) 


1287r 2 a 2 
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The terms proportional A m 2 ip 2 and A 2 /4 ip 4 in (g9|) are cancelled by the mass and coupling constant counter terms. 
The divergent term which depends on m 4 determines the cosmological constant counter term, that is 


6A = 


m 


64-7T 2 1 e 


—h In 


47r/i 2 


(90) 


The remaining terms in (89) combine with the corresponding expressions of the following two dimensionally regularized 
integrals 


-5(6£-l)(f +H 2 ) j 


1 


1 


d 3 k 


(2tt) 3 L2S^ fc0 4 Cl 3 


-V(r) 


kO 


trr ^( R , „2^ M 2 (t) (2 47T/i 2 

= (6«-i) (e + H jlSv{7 + ln Myor 7 

+ 3^< 6 f- 1 >(f + H2 ) M2(0) - 


(91) 


and 




a J (27r) 3 4f2| 


ko 




1 


32tt 2 


(6£ — 1)— [(£ — |) R'a 2 + A (pip' — \aH<p 2 ] | - + In 


47r/i 2 

M 2 (0) 


(92) 


The term A(£ — 1 /6)Rtp 2 in ( |89|) is cancelled by the same term with opposite sign in (|9l|). The H 2 ip 2 -terms in 
and ( p2| ) are absorbed by the counter term proportional to 8 £, and the divergence proportional to ipip 1 is compensated 
by this counter term as well. The remaining (^-independent but still time-dependent divergent terms are absorbed 
into the counter terms 5a H t t and 5Z G t t■ We choose 


32t r 2 


5a = — 


5Z = — 


ATTfJ 2 


- + In 
__ e nv 

(^- s)w 2 f 2 


47 TfJ 2 

m —-- 7 


(93) 

(94) 


That is, all divergent integrals appearing in the unrenormalized divergent fluctuation integral of the energy are removed 
by the corresponding counter terms and we finally arrive at the renormalized expression for the energy 


1 


1 


A +AA 


£ ™"=2^*’' 2 + 2^ ( "* 2 + Am > 2 + 4!a* 


- 6 (£ - i + AC) 


H 


H 


2 ^ ^ 


+AA + A aH tt + AZG tt + |^0(6£ - 4) (f + H 2 


M 4 ( 0) M 2 (0)M 2 (r) 


.1 f d 3 k ( V(t) / 
a 2 J (27t) 3 | 4flfeoa 2 V 


2Reftf+ |fc™| 2 ) + 3n -^K 


l^l 2 ) 


1 


1287r 2 a 4 
( 2 ) ,2 


327r 2 a 4 


1 


2f2fcoa 2 

1 


Re 


K^’) + Tii4e c2(V "' r) + T«ib 52(F ' T) 




i? 


90 2 


C(W,r) + 2Re/i fc ®+ |4 1} | 2 


+^- (6 £- 1) - 
4* & 


1 C(V", r') + -L- cos 2fWW(0) 


2^1-0 




+2Re^ 2) ' +2Re4 1) *^ 1)/ 


(95) 


The finite corrections are 
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Aa = 


(96) 


(f-E) 2 h i. ™ 2 


32t r 2 M 2 (0) ’ 


4 2 

A ~ m , m 
AA = -In- 


AZ = 


64tt 2 M 2 (0) ’ 


2 9 

, m~ m 

— -In- 


167T 2 


M 2 ( 0) 


(97) 

(98) 


Next we have to consider the renormalization of the trace of the energy momentum tensor. We introduce the 
available counter terms into the unrenormalized expression for T so that it reads now 


T£ = -[i-6(£ + ao] 


. _/2 tt2 

V , H ..2 


R 


+ + 


6(S + 50 


-w 


+2(?n 2 + 5m 2 ) + 45A + S&m + 5ZG() 

a 2 6a 4 ^ ^ 


1 r d 3 k j 1 

a* J (2^)3 | a 5 


(1 - 60 


1 


1 


£7 


2f2fco 

d 


1 


l^ (T) | 2 + | (~2*Imh' fe (T) - 2iIm4 T) *^ (T) ) 

y(r)(l + 2Re/i fc (T) + |4 1) | 2 ) 
l + 2Re/ijfc (T) + |/if } | 2 ) 


2ft 


/c0 


"* 2 + S?*’" - (1 “ 60 ( ff2 “ 


R 


x (l + 2Reh fe (1) + |/ii 1) | 2 ) | ■ 


(99) 


We can split the trace of the stress tensor into a divergent and into a convergent part. The divergent part reads 


T a* _ 

A 4 div 


d 3 k 1 


(2tt)3 


(1 - 60 


8(fifco)' 


■ V"{t) - (1 - 60 


H 1 
a 4(f2fco)’ 


:V\t) 


+ 


m- + —.tp 2 - (1 - 60(R 2 - gR) 


2a 2 


1 


1 


2£O 0 4(f2feo)' 


W(7 


( 100 ) 


The first derivative of the potential cancels with the corresponding term from ( |50| ) in dl00|) . The remaining divergent 
parts are 


1-60 

n™ n 2 / (27r) 3 i 1 ^ 


8(Ofco) 3 


\RM 2 {t) - 2 H 2 M 2 {t) 


H 


+2aH (£ - i) R' - 2A- w ' + (£ - ±) R 


+^\{W + v?' 2 ) - A|R^ 2 + \H 2 ip 2 


” 2+ ^ 2 -< i - 6 «( fl2 -! 


L 2S7fc 0 4(f0o)' 


:V(t) 


( 101 ) 


After dimensional regularization again the counter terms absorb all divergent terms in the fluctuation integral of Tjf. 
The renormalized trace of the stress tensor takes the final form 


TH = -[ l-6(£ + A0] 
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+2 (m 2 + Am 2 )^ + A ^ f''V + 4AA + A.aH£ + A ZG* 


M 2 { 0) 

167r 2 a 2 

1 f d 3 k I 1 


o' 6a 4 5 
2 A 
2 a 2 


^■ + ^ 2 -(l-6?) U 2 - 7 


a 2 J (27r) 3 I a' 


(1 - 60 
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2 Q fco K (1) | 2 + | (—2ilm /i^ 2) - 2*Im4 1) *ftfc 1) 
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V(t) (2R eh k ™ + |4 1} | 2 ) + ^2-F'(0)sin2f) fc0 T 


fc k 

1 

2 

kO 


— -V"(0)cos2O fc0 r- —-C{V"',T) 


+ 2fi fc0 a (1 6C) 


—-C(U",r)+2Re/^ 


/( 2 ) 


90 2 

Zi */c0 


cos2D fc0 TU'(0) + 2Re/i , fc (1) */ijj. 1) 


212 


A:0 




1 


20 1 2 
2:1 ‘fcO 


C(F , ,r)+2Re/i fc (2) + |4 1) | 2 


( 102 ) 


VIII. CONCLUSIONS 

We have presented here the renormalized equations of motion for a scalar field in a conformally flat FRW universe 
including the quantum back reaction in one-loop approximation. We have used dimensional regularization and an 
MS renormalization. However, within this formalism other renormalization conditions can be employed, as well as 
other regularizations. The formalism is fully covariant and the counter terms can be chosen independent of the initial 
conditions. In contrast to the adiabatic regularization it is not based on the WKB expansion and can therefore be 
generalized to coupled systems ||29| ] . Furthermore, we do not have to perform delicate subtractions in the divergent 
integrals since in our formulation they are finite from the outset. The method can be adapted easily to finite 
temperature computations as well, since no new ultraviolet divergencies are introduced. 

The intention of this work is of course to supply at the same time a useful scheme for numerical computations. We 
have shown previously ||28| , |29| that the numerical implementation is straightforward. A realistic application to one of 
the inflationary scenarios needs, however, a judicious choice of initial conditions and of the parameters of the theory, 
as well as extensive numerical experiments. While such computations are in progress we prefer to present here the 
formalism as such; we think that is by itself of general theoretical as well as of practical interest. 
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